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Abstract 

Concepts from Ergodic Theory are used to describe the existence of special non- 
transitive maps in attractors of phase synchronous chaotic oscillators. In particular, 
it is shown that for a class of phase-coherent oscillators, these special maps imply 
phase synchronization. We illustrate these ideas in the sinusoidally forced Chua's 
circuit and two coupled Rossler oscillators. Furthermore, these results are extended 
to other coupled chaotic systems. In addition, a phase for a chaotic attractor is de- 
fined from the tangent vector of the flow. Finally, it is discussed how these maps can 
be used to a real-time detection of phase synchronization in experimental systems. 
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1 Introduction 



Coupled chaotic systems are recently calling much attention due to the veri- 
fication that they may be useful to the understanding of natural systems in a 
variety of fields as in ecology p], in neuroscience [2|3] , in economy [I], and in 
lasers [5|6j . It has been verified that despite of the higher dimensionality of a 
coupled chaotic system, the coupling among the elements might make them 
to synchronize [7f5] . reducing the dynamics of the system to a few degrees of 
freedom. 
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In this work, we focus our attention in the phenomenon of Phase Synchroniza- 
tion (PS), which describes the appearance of a phase synchronous behavior 
between two interacting chaotic systems [5J, i.e., given two chaotic systems, 
their phase difference remains bounded, despite of the fact that their ampli- 
tudes may be uncorrelated. This phenomenon is particularly interesting since 
it can arise from a very small coupling strength. Its presence was reported 
in a variety of experimental systems. It was experimentally demonstrated in 
electronic circuits [ID] , and latter in electrochemical oscillators [IT]. It was 
found in plasma [12], in the Chua's circuit JT3], and there were also found 
evidences of phase synchronization in communication processes in the Human 
brain [HfTS J . 

To detect PS in a real-time experiment, one has to measure the phase of the 
chaotic trajectory [IB]. However, the phase is not always an easily accessible 
information. To overcome this difficulty, it is important to understand funda- 
mental properties of phase synchronous systems, that could be experimentally 
easily verified. For chaotic systems that are phase synchronized with a periodic 
forcing [IT], it was reported that a stroboscopic map of the trajectory was a 
subset of this attractor and occupies only partially the region delimitated by 
a projection of the attractor. This property was used to detect in a real-time 
experiment phase synchronization between the Chua's circuit and a sinusoidal 
forcing [13] . 

This approach of detecting phase synchronization through the stroboscopic 
map can be extended for coupled chaotic oscillators, in a formal way. The 
stroboscopic map is generalized to the Conditional Poincare Map. Given 
two oscillators, at least one being chaotic, the conditional Poincare map is con- 
structed by collecting points in one oscillator at the moment at which some 
event happens in the other one. If the set of discrete points generated by this 
conditional map does not visit any arbitrary region of a especial projection 
of the chaotic attractor, we call this set a P-set. This property of the con- 
ditional Poincare map is called non-transitivity [18], i.e., an initial condition 
under the conditional Poincare map does not visit everywhere in a subspace 
of the attractor. Alike the stroboscopic maps of oscillators that are in phase 
synchrony with a forcing, the conditional Poincare maps of coupled chaotic 
oscillators, in PS, also only partially occupy a projection of the attractor. 

In this work, we show how the conditional Poincare map can be used to 
detect PS, without actually having to measure the phase. For phase-coherent 
oscillators, a special type of P-set, that we call PS-set (Phase Synchronization 
set), exists. Conversely, its existence also implies PS. We illustrate our findings 
and ideas with numerical and experimental analyzes in the forced Chua's 
circuit, and the coupled Rosller oscillator [19J. 

Further, we extend these results to non-phase coherent attractors. Finally, 
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we also introduce a phase of a chaotic trajectory to be a quantity related to 
the amount of rotation of the tangent vector. This definition can be used to 
chaotic attractors, independently whether they have phase-coherent or non 
phase-coherent dynamics. 

This work is organized as follows: in Sec. (2J we define a way to measure the 
phase of a chaotic flow, and discuss the relation between the average returning 
time and the average angular frequency. In Sec. [3], we discuss the conditions 
for PS states and, in Sec. HI we describe the phenomenon of PS in the forced 
Chua's Circuit. We introduce the notion of a conditional Poincare map in Sec. 
and the P-sets (as well as the PS-sets) in Sec. El In Sec. [TJ we show how 
PS can be found by the detection of these sets in the forced Chua's Circuit 
and in Sec. [H for the coupled Rossler oscillator. Further, in Sec. M> we discuss 
the extension of these ideas to a class of non-coherent oscillators. In Sec. [101 
we make some remarks and the conclusions of this work. In Appendix A, 
we formally introduce the conditional Poincare map and the P-set, and in 
Appendix B, we show that for coherent dynamics the PS-sets exist if, and 
only if, there is PS. In other words, PS implies PS-sets and vice-versa. 



2 Phase, frequency and average returning time of a chaotic attrac- 
tor 



The phase of a chaotic attractor in a projection j (a subspace) is defined to 
be the amount of rotation of the tangent vector in this projection, and can be 
represented by an integral function of the type 



/' 


d6(t') 


/o 


dt' 



dt 



with d9{t) being an infinitesimal displacement of the tangent vector of the 
flow, from time t to time t + dt, and dt — > 0. Note that in Eq. flTJ), we are 
measuring the amount of rotation, per unit time, of a projection of the tangent 
vector of the flow, on the same subspace j where the phase is defined. We call 
this subspace Vj. The attractor X, projected on the subspace Vj is regarded 
as Xj. The instantaneous angular frequency of the trajectory in Xj, named Wj 



is given by So, from Eq. ([T]), Wj 

(Wj) is 

m - 



dt 



and, the average angular frequency 



dt ' 



(2) 



(— ) represents the average. Equation 

0jW 



can be put into the form (Wj) 



Whenever a Poincare section can be defined, the average period of the chaotic 
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attractor on the subspace Vj is calculated by 



(Tj) = (3) 

where Arj = r- — Tj^ 1 , and rj represents the time at which the trajectory in 
the subspace Vj makes the i-th crossing with this Poincare section. 

We introduce (A0j) to be the average displacement of the phase for a typical 
period as 

(A*, ( ) = 4f>, (4) 

with cj)j(N) being the phase associated to the subspace Vj at the moment that 
the N-th crossing between the trajectory and the Poincare section happens. 



Thus, we can put Eq. (jSJ) as 



For the forced Chua's circuit, the subspace V\ is defined by a suitable pro- 
jection of the circuit variable. We have that (A0x) = 2n. So, Eq. (jSJ) can be 
written as (W\) = 7§ry. For the coupled Rossler oscillator, the quantities in Eq. 
(jSJ) can be calculated in two subspaces, one associated to the variables of one 
Rossler, the subspace Vi, and the subspace V2, associated to the variables of 
the other Rossler system. As shown in |2Tj, might slightly differs from 

2vr, and thus, (Wj) = 

So, Eqs. (j2J) and fl5j) relate the average period, the average angular frequency, 
and the phase of a chaotic trajectory. This shows that the average period 
(recurrence) and the average angular frequency are intimately connected in 
phase coherent chaotic systems, and both these quantities can be calculated 
from the phase. 



3 Phase Synchronization 



Having defined phase, PS exists whenever the following condition is satisfied 



\Mt)-rMt)\<(^i)- (6) 

The minimal bound for the phase difference, (A0i), in terms of the phase 
as defined by Eq. dSJ), was theoretically estimated in Ref. [21]. Equation (J6]) 
means that the phase difference between the two coupled systems is always 
bounded, and r is a rational constant [2"2"] . 
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Also, 

(W 1 )-r(W 2 ) = 0. (7) 
In this work, we will consider cases for r=l. Otherwise, a simple change of 
variables could eliminate this constant from Eq. ATI) . For the forced Chua's 
circuit (W2) = u), with uj = 2irf, representing the angular frequency of the 
forcing. There is PS, if (Wi)=u. Therefore, (Xi)=l//. For the coupled Rossler, 
if PS exists, (Wi) = (W 2 ), (R 1 ) = (R 2 ), and (A0i) = (A0 2 ), and therefore, we 
could have in the right term of Eq. (jfJJ) (A02), instead of (A0i). 



4 The sinusoidally forced Chua's circuit 



The circuit is represented by: 



Ci-^ = <7(X 2 -*i)-i*L (8) 
jy 

C"=g(X 1 -X 2 )+X 3 (9) 
at 

L-^-i = -X 2 - V sin (oot) (10) 

where X 1; X 2 , and X 3 represent, respectively, the tension across two capacitors 
and the current through the inductor (See [13J for more details), uj and V are 
the angular frequency and the amplitude of the forcing, respectively. The 
piecewise- linear function, i NL , is given by: 



i NL = m X l + 0.5(mi - m )[|Xi + B p \ - \X X - B p \] (11) 
where we have chosen the parameters C\=0.1, g=0.574, C 2 =l, L=l/6, m =- 
0.5, and m^-0.8, such that we obtain a Rossler-type attractor, for V=0. 

To calculate the phase of the chaotic trajectory, we first define the subspace V\ 
to be given by the pair of variables (X%, X 2 ), and then, we use Eq. JTJ. In Fig.[[J 
we show the difference between the phase of the chaotic circuit [as calculated 
by Eq. ([rj] and the phase of the forcing, out. In (a), the phase difference 
is bounded and the average period of the chaotic attractor, defined as the 
average recurrence time of trajectories that cross the section X 2 = 0, is equal 
to (Ti)=3. 57015, which is equal to 1/f, since /=0.2801. The average angular 
frequency can be calculated using Eq. (jSJ), which gives us (Wi)= 1.75992. Or, 
from Eq. (T5]), which gives us (Wi)—1. 75992. Note that the average growing 
of the phase [calculated by Eq. ((T])] for a typical average period is 6.28318 . . . 
which is 2tt. In Fig. [T](b), we have that (Ti)=3. 57006 which is different from 
1/f, since /=0.279. So, in (b) there is no PS, and consequently Inequality (jSj) 
is not satisfied. 
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time 



Fig. 1. (a) Phase difference is always smaller than 2ir, so PS exists between the 
circuit and the forcing, for /=0.2801 and V=0.0015. (b) PS is not present and the 
phase difference grows bigger than 2tt. 

In Ref. [1 3 j , we detected PS experimentally in the forced Chua's circuit, 
whenever stroboscopic maps could be constructed for a time interval equal 
to Atx = j, such that this map, projected into the same subspace considered 
to calculate the phase, does not occupy the region occupied by the attractor 
projected on the same subspace. 

To understand this technique, we assume that <pi(t), the phase of the chaotic 
trajectory, is the angle (on the lift) described by the vector position of this 
trajectory (Rossler-Like attractor), and 2 ()f:)=co , t the phase of the forcing. If 
Eq. ([6]) is satisfied at any time, then it is satisfied at multiples of the period 
of the forcing r\ = j. So, we get \<f>i(rl) — 2ni\ < 2n, which means that a 
stroboscopic map has to be concentrated in an angular section smaller than 2n. 
The stroboscopic map, which is already a subset of the chaotic flow, projected 
into the same subspace considered to calculate the phase, does not occupy 
the whole region occupied by the attractor projected in this same subspace. 
Another property of the stroboscopic map is that points in it are mapped into 
it by looking at the trajectory after a time interval given by Ati, so, it is a 
subset that is recurrent to itself. 

Using this technique, and for the same parameters as |13] . we show in Fig.[2^a) 
the experimental synchronization region for the forced Chua's circuit in the 
parameter space V xf. The triangular shaped region represents parameters for 
which the stroboscopic map has the points concentrated in an angular section 
smaller than 2n. The bump at the bottom right side of the PS region is due 
to non- synchronous states that present a long bounded phase difference. By 
a typical time interval within which the experiment is realized, which was of 
the order of 40,000 cycles, the system seemed to be phase synchronized, i.e., 
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localized stroboscopic maps were found. In fact, we have detected these maps 
even for observation times corresponding to 150,000 cycles, in the region of 
the bump. 

In a short, the bump region is an extended structure in the circuit parameter 
space, that presents intermittent behavior in the phase difference [24J, but 
with a long laminar regime, even for parameters far away from the border 
between the PS and the non PS region. This intermittency differs from the 
usual one, observed in the transition to PS, in the fact that this latter happens 
very close to the border between the PS and the non PS region. The reason 
for this intermittency is due to the presence of a periodic window, close to the 
region of the bump. 

Simulation is shown in Fig. 12(b), where black points represent perturbing pa- 
rameters for which Eq. ([6]) is satisfied, with <fti defined in Eq. ([1]). One sees 
that the PS region resembles a triangle. The triangular shaped region, de- 
noted by the light gray dashed line, represents the region where the system is 
not phase synchronized, but the phase difference remains bounded for a long 
time interval, that might be longer than 100,000 cycles of the systems. So, we 
reproduce numerically the same atypical intermittency observed experimen- 
tally, i.e., long laminar regime in the phase difference, for parameter regions 
away from the border between PS and non PS states. This happens associated 
with a periodic window, as the one shown in Fig. [21(b). 




Fig. 2. (a) Experimental PS parameter space, (b) Simulated PS parameter space. 
Black points represent parameters for which Eq. ([6]) is satisfied for 120,000 crossings 
of the trajectory at X2=0. In both figures, the horizontal axis represents the forcing 
frequency / and the vertical axis its amplitude V. Variables in (b) are dimensionless 
and /o is the main frequency of the non- forced circuit. 

The shape of the synchronization region in Fig. [21(b) is equivalent to the region 
in the experiment, constructed by detecting the stroboscopic maps contained 
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within a small angular section. This proposes an equivalence between the 
existence of a recurrent subset and the verification of Eq. . Inside the syn- 
chronization region, a stroboscopic map appears like in Fig. E](a), where the 
light gray points represent the attractor, and the dark filled circles, the map. 
Outside of the PS region, there are parameter sets for which the stroboscopic 
maps do not occupy the region occupied by the attractor, at the projection 
in which the phase is calculated. As one sees in Fig. [3]^b) [for the parameters 
represented by the plus symbol in Fig. Efb)] there is a region of the projected 
attractor (pointed by the arrow) for which the stroboscopic map never visits. 

The difference between the stroboscopic map that appears while there is PS 
in (a), and the stroboscopic map that appears while there is not PS in (b), is 
that points of the stroboscopic map in (a) are all concentrated in an angular 
section of the attractor. As it will be further classified, the stroboscopic map 
in (a) is a PS-set and in (b) is a P-set. 




Fig. 3. A projection of the attractor in light gray, and its stroboscopic map in 
dark gray. The parameters are ,f =0.2801 and 1^=0.0015, in (a), and /=0. 28063 and 
^=0.0016, in (b). 

As a way to better characterize the PS phenomena in the Chua's circuit, we 
calculate the Lyapunov exponents. As usually expected, the transition to PS 
is associated to one of the exponents becoming smaller than zero. Since there 
is already one exponent that is smaller than zero, PS induces the creation of 
a second stable direction in the Chua's circuit. Previous to the transition to 
phase synchronization, this exponent was zero. In Fig. m we show the exponent 
(and the error bars) associated to PS in black and the quantity (Wi) — uj (in 
gray, for a fixed amplitude of \^=0.0015) with respect to the frequency. In the 
region that the exponent becomes smaller than zero, (Wi)=u, satisfying Eq. 

©• 
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Fig. 4. The Lyapunov exponent (and the error bars) associated to PS in black, and 
the quantity (Wi) — u in gray, for a fixed amplitude of ^=0.0015 and a varying 
frequency. 

5 Conditional Poincare Map 

The finding of maps of the attractor that appear as localized structures implies 
PS. The conditional Poincare map introduced in this chapter as a generaliza- 
tion of the stroboscopic map, is an efficient way of revealing the existence of 
such special mappings. 

The stroboscopic map technique defined for periodically driven chaotic sys- 
tems was explained in the previous sections. To extent the idea of stroboscopic 
map in coupled chaotic oscillators we came up with the conditional Poincare 
map, which is a map of the flow, constructed by observing it for specific times 
at which events occur in the subsystems X y An event is considered to happen 
when the trajectory of one subsystem crosses a Poincare section. So, given 
two oscillators S\ and S2 (at least one being chaotic), with trajectories in the 
subspaces where the phase is defined, the conditional Poincare map of S\ is 
the trajectory position at the moment that a series of equal events happens, in 
S%. Analogously, the conditional Poincare map of S 2 is the trajectory position 
at the moment that a series of equal events happens, in Si. In the case of 
periodically forced chaotic systems, an event may be defined to happen when- 
ever the forcing reaches a specific value, and the conditional Poincare map is 
the usual stroboscopic map, because the time interval between two events is 
constant. In coupled chaotic systems, the time interval between two successive 
events is no longer constant. 

We define a time series of events r], by the following rule: 
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• t\ represents the time at which the i-th crossing of the trajectory of £2 
occurs in a Poincare plane. 

• t\ represents the time at which the i-th crossing of the trajectory of Si 
occurs in a Poincare plane. 

The discrete set of points observed at the times rj is called set V. This set, 
projected at the subspaces Vj (where the phase is calculated), is named Vj. 
The conditional Poincare map is represented by T t j . So, we say that Vj is the 
set of points generated by T T i . 

The next step is to define when Vj can be regarded as either a P-set or a 
PS-set, this last set implying phase synchronization. 



6 Sets generated by the Conditional Poincare Map 

The T>j set is a P-set, if it does not completely fulfill the projection Xj of the 
attractor. In other words, a discrete set Vj is considered to be a P-set, if for 
balls of radius 5, centered in all points of the attractor projection Xj, one does 
not find points of Vj inside of all these balls. If Vj completely fulfill Xj, we 
say that these two sets are equivalent (and we represent this by the symbol 
=). More details, see Appendix A. 

So, a P-set exists, if the conditional Poincare map is not transitive on Xj |18j . 
That is, the flow, observed by the times for which the conditional Poincare 
map is defined, does not visit arbitrary regions of Xj. Note that however, the 
attractor is chaotic, and therefore, the chaotic set is always transitive through 
the flow. So, given a set of initial conditions, its evolution through the flow 
eventually reaches arbitrary open subsets of the original chaotic attractor. 

We can classify three relevant types of sets generated by the conditional 
Poincare map 

type-a Vj is equivalent to Xj (Vj = Xj). The conditional Poincare map T T j 
is transitive in Xj. 

P-set Vj is NOT equivalent to Xj (Vj ^ Xj). The conditional Poincare map 

is NOT transitive in Xj. 
PS-set Vj is P-set, with the additional condition that it is localized in the 

vicinity of the Poincare section chosen to define the events. 

In the following, we comment each case: 
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6.1 type- a sets 



These sets appear whenever there is no PS. If two non-identical coupled chaotic 
systems (topologically similar) are not phase synchronized, the chaotic tra- 
jectories do not make correlated events in both subspaces V\ and P 2 - As a 
consequence while the trajectory is positioned at the specified Poincare plane, 
at the subspace Vi, the trajectory in the subspace V2 is everywhere in this 
subspace, making the set Vj to be equivalent to Xj. 

An interesting illustration is the case of two uncoupled equal chaotic systems, 
but with different initial conditions. As we construct the conditional Poincare 
maps, they will be a type-a set, since also the distance between the trajectories 
in the two oscillators are sensitive to the initial conditions, and will diverge 
exponentially. 



6.2 P-sets 

These sets constitutes a necessary, but not sufficient, condition to describe PS. 

In periodically forced chaotic systems, they might appear when there is not 
PS. As we already mentioned, the points in these sets are not localized in 
special spots of the attractor projection. As a consequence, the domain of 
the absolute difference between the time at which the same number of events 
happen in both oscillators has a broad character. 



6.3 PS-sets 

PS-sets imply PS and vice-versa. They exist, if and only if, there is phase 
synchronization, as shown in Appendix B, and illustrated with the examples 
in Sec.Hl and throughout this work. Another important point is that the PS-set 
provides a real-time detection that can be easily experimentally implemented 
(Sec. HJ), and easily constructed from a data set. 

PS-set implies PS because the difference between the time at which the iV-th 
event happens in both oscillators is small, which means that the time difference 
— 7vf I is smaller than a finite constant value. As a consequence, the points 
in the conditional Poincare map of one oscillator are confined around the 
Poincare section chosen to define the events. Therefore, the detection of a PS- 
set can be done by observing this characteristic of the conditional Poincare 
map. 
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For Rossler-like oscillators, in which the trajectory spirals around an equilib- 
rium point, the PS-set is confined within an angular section. 



6.4 Length of a PS-set 



Having found the PS-set, we can study properties of these sets that give us 
the level of organization and coherence of the oscillators 

A PS-set, T>j, is said to have length 1, if the set is constructed by the time 
series of the same events. Defining the event to be given by the crossing of the 
trajectory to a Poincare plane, the corresponding time series of events Mj(l) 
is given by rj, t] +1 , t] +2 , t] +3 , . . .. For this PS-set, points in Vj are mapped in 
Vj, after one application of the conditional Poincare map. 

A PS-set is said to have length 2, if it is obtained by a time series of two 
different events. So, a PS-set can be constructed from more complex series of 
events. We construct a length-2 basic set using a time series of events Mj(2) 

given by rj, - J — ^ — , rj + , J — ^ — , . . . [25]. As an example, for the perturbed 
Chua's circuit, the length-2 basic set is constructed by a stroboscopic map that 
collects points every half period of the forcing. A length-2 basic set is assumed 
to be composed by two other subsets, named minimal sets, the subsets V® and 
V). 

They have the property that if a point x is such that x G £>°, this point, 
iterated by the conditional Poincare map, goes to the minimal set Pj, and if 
xo is such that x 6 Vj, this point, iterated by the conditional Poincare map, 
goes to the minimal set Z>j. Thus, points in "Dj are mapped to itself after 2 
applications of the conditional Poincare map. The minimal set £>j is said to 
be disjoint to the set if they do not intersect, i.e., X>j fl D 2 - — 0. For some 
systems that present a strong phase-coherent state, as the ones here studied, 
in which the instantaneous trajectory velocity does not differ too much from 
the average velocity on typical orbits, it is possible to find length-2 basic set, 
with disjoint minimal sets, when PS is present. 

For a general case, we do not expect to find a length-2 PS-set with disjoint 
minimal sets. As an example, one can think of a spiking-firing oscillator, phase 
synchronized with a periodic forcing. Due to the fact that the spiking-firing 
dynamics has a fast and a slow mode, the conditional Poincare maps might 
overlap. 
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6.5 Sets diagram 



Here, we explain through a diagram the possible emerging sets from the con- 
ditional Poincare map. 




type - a P - set 

PS < > PS-set 

Starting from the chaotic attractor X, the set T> is constructed from the condi- 
tional Poincare map, represented by T, we project V and X into the subspace 
Vj, obtaining the sets Vj and Xj, respectively. 

We classify the set T>j into type-a set or P-set, by checking whether the con- 
ditional Poincare map is transitive in Xj, i.e., by verifying the equivalence 
between the sets Xj and T>j. Then, if the P-set is localized in the vicinity 
of the Poincare section where the events occur, the P-set is a PS-set, which 
means that PS is present. 



7 PS-sets in the Chua's circuit 

For applying our formalism to the periodically forced Chua's circuit, the event 
times are r\ = ir, with r representing the forcing period. The time series for 
the length-2 basic set is given by l/2r, r, 3/2r, 2r, 5/2r, 3r, . . .. Everywhere 
inside the PS region, we find length- 1 [as an example, see Fig.[3](a)] and length- 
2 [as an example, see Fig. 0(a)] PS-sets, this latter with disjoint minimal 
sets. In Fig. 0(a), the application of the conditional Poincare map in the 
minimal set V>\ leads to the minimal set T>\. Both sets form a length-2 PS- 
set. Note that both T>\ and V\ can be regarded as a length- 1 PS-set. To our 
numerical precision, we have checked that there are not basic sets beneath the 
Synchronization region tip. Which means that type-a set is present for very 
small but finite amplitude forcing. Outside the PS region, there is a P-set, i.e., 
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a non-transitive conditional Poincare map on the chaotic attractor projection. 
So, T>j ^ Xj. In this case, there is not PS. More examples of length-Q PS-sets 
in the phase synchronous forced Chua's circuit can be seen in |26| . 




X 2 



Fig. 5. (a) Length-2 PS-set with disjoint minimal sets, in the forced Chua's circuit, 
for the parameters / = 0.2801 and V = 0.001. The PS-set is constructed by the time 
series Mi (2). (b) Length-2 PS-set with disjoint minimal sets, in the coupled Rossler 
oscillators, for the parameters e=0.01 and 5a = 0.001. The PS-set is constructed by 
the time series M^{2). 

8 PS-sets in the coupled Rossler system 

We can use the formalism of the conditional Poincare map to study the appear- 
ance of PS in coupled chaotic systems, as the two coupled Rossler oscillators 
given by: 

a?i,a = -«i,22/i - Z\ + e{x 2> i - x 1>2 ) (12) 
2/1,2 = "1,2^1 + 0.15j/i (13) 
ii, 2 = 0.2 + Zl (xx - 10) (14) 

with ai — 1; and a 2 = a± + 5a. The index denotes systems 1 and 2. The 
subspaces Vj are defined by Vj—(xj,yj). In a coupled chaotic system, r\ (j % 2 ) 
does not increase uniformly, but it is given by the time the trajectory crosses 
the Poincare plane j/2=0 (2/1=0) . For these times, and using the time series 
Mj(2), we construct the minimal sets T>\ and T>\ (T>\ and T>\). The parameters 
are e=0.01 and 5a = 0.001. 

In Fig. [5(b), we show a length-2 basic set with disjoint minimal sets. The 
application of the conditional Poincare map in the minimal set T>\ leads to the 
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minimal set T>\ y and vice- versa. These two sets form together a length-2 PS-set, 
but each one separated can be regarded as a length-1 PS-set. The characteristic 
of these PS-sets is that they appear as localized structures around the Poincare 
section chosen to define the events. As one can see, the set T>\ is localized in 
the neighborhood of the line yi = 0, where the Poincare section is chosen. 

In fact, this length-2 PS-set with disjoint minimal sets (as well as a length-1 
PS-set) is found everywhere in the PS region, as shown in Fig. O In it, filled 
squares represent parameters for which these special PS-sets are found, and 
empty circles parameters for which PS exists. 

A PS-set of length-Q is detected using in Eq. flA~2|) V® \ from which we can 
check whether occupies (type-a discrete set V) the whole space occupied 
by X 1 . The set Bi(x) in Eq. f1A.2h is constructed assuming squares of size 
I = 1.5, in points of the set V^ 1 - In Fig. El we show a case for Q=2. 



0.0015 - 



OJ 0.001 - 





Fig. 6. Empty circles show parameters for which PS exists, detected by Eq. ([6]), and 
filled squares represent parameters for which a length-1 PS-set appears simultane- 
ously with a length-2 PS-set with disjoint minimal sets. Horizontal axis represents 
parameter mismatch and vertical axis, the coupling amplitude. 



In Fig. [3(a), we show the attractor in the subspace V\, the subset X\ in gray, 
and the discrete set T>i in dark empty circles . In (b), we show a magnification 
of the box in (a). Note that neighborhoods of arbitrary points in the trajectory 
of Xi (gray) always contain a point of the discrete set T>i (dark empty circles). 
So, the set X x is equivalent to the set T>i and, therefore, T>\ is not a PS-set, 
and therefore, it does not exist PS. Differently from the Chua's circuit, the 
Rossler coupled system presents no P-sets for parameters outside of the PS 
region. 
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Fig. 7. We show a situation where PS is not found for e=0.01 and 5a = 0.00001. 
In (a), we show the attractor in the subspace V\, the subset X\ in gray, and the 
discrete set T>\ in dark empty circles. In this figure, T>\ is a type-a set, since X\ visits 
every neighborhood of points in T>\ (the conditional Poincare map is transitive in 
X\). In (b), we show a magnification of the box in (a). 

9 Extension to other coupled chaotic systems 

The approach presented in this paper can be extended for non-coherent at- 
tractors. As noticed in Ref. [20], attractors that present non-coherent phase 
motion in the phase space may present a coherent motion in the space of the 
velocities, i.e. (x, y), which is the case of the funnel attractor [20]. In this case, 
the extension of our approach to the non-coherent phase motion is straight- 
forward. Instead of defining a conditional Poincare map in the phase space, 
we analyze the dynamics in the velocity space, in which the phase is coherent, 
and therefore, all the ideas introduced herein can be applied in this new space. 

To some chaotic attractors, it might not be possible to define a Poincare section 
to construct the conditional Poincare map. However, one can still construct 
these conditional maps by denning different types of events. For example, in 
coupled neurons, this event can be chosen as being a beginning or the ending 
of the bursts, which can be well denned by the crossing of the trajectory with 
some given threshold. 
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10 Conclusion and Remarks 



A chaotic set is always transitive through the flow. So, given a set of initial 
conditions, its evolution through the flow eventually reaches arbitrary open 
subsets of the original chaotic attractor. However, a stroboscopic map of the 
flow, whose generalization is here called conditional Poincare map, might not 
possess the transitive property. That is, given a set of initial conditions, its 
evolution through the conditional Poincare map might not reach arbitrary 
open subsets of the chaotic attractor. 

The introduction of the term "conditional" in the map nomenclature comes 
from the unconventional and non rigidly way we adapt the established def- 
inition of a stroboscopic Poincare map. For coupled chaotic oscillators, this 
conditional map is constructed based on events, conveniently chosen at the 
same subspaces where the phase of a chaotic system is defined. In addition, 
the application of this map through the flow, which results in a discrete set, is 
inspected not in the whole phase space, but also in the same subspaces where 
the phase is defined. 

If phase synchronization exists the conditional map generates special discrete 
sets, named PS-sets. The contrary is also true, i.e., a PS-set implies phase 
synchronization. This was illustrated in the periodically forced Chua's circuit, 
in the coupled Rossler oscillator, and in other more general (topologically 
equivalent) coupled chaotic systems. 

The ideas introduced here provide an efficient way of detecting phase syn- 
chronization, without having actually to measure the phase of the chaotic 
trajectory. Indeed, this detection can be done in experiments in a real-time, 
as it was done here in the perturbed Chua's circuit. 

It is worthy to say that the PS-sets are robust under small additive noise that 
could corrupt the data in an experiment. This is so, because a small additive 
noise does not interfere much with the time that the trajectory crosses the 
Poincare section, but just deviates the timing of the crossing. Then, the PS- 
sets still remains. This robustness against the noise is an important property 
in order to apply these ideas in experiments, as done in this work. 

We have also introduced the phase as a quantity that measures the velocity 
of rotation of a projection of the tangent vector along the trajectory. This 
definition can be used to arbitrary flows, independently whether or not they 
present a coherent or non-coherent phase dynamics. 

Finally, our formalism of the conditional Poincare map can be used in coupled 
maps (or perturbed) to detect synchronous behavior (not phase synchroniza- 
tion) between the systems (or between the forcing), for the case where one 
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does not find full synchronization between the maps. One particular example 
where that happens is in the periodically forced Logistic equation [27] or for 
a system of coupled Logistic maps [28j, where one can find a finite number of 
synchronous chaotic subsets, the basic sets. 

Acknowledgment Research partially financed by FAPESP, Alexander Von 
Humboldt Foundation, and CNPq. 



A The conditional Poincare map and the Vj sets 

Next, we present some formalism using elements of Ergodic Theory in order 
to introduce the conditional Poincare map and the Vj sets. 

Given a flow F t , we call X the chaotic attractor. A subset of the attractor is 
its projection on the subspace Vj, that we call Xj. We assume X e R d . 

The notion of stroboscopic map of a flow can be generalized to any tempo- 
ral translation on the trajectory. We define the temporal translation to be a 
transformation represented by T T i , called conditional Poincare map. The ini- 
tial condition x(t) is iterated under the temporal transformation T t j , to the 
point x(t + Tj). Applying the transformation T t j in a typical trajectory, for 
the time sequence rj, Tj +1 , rj +2 , give us the points x(rj), x(rj +1 ), and x(rj +2 ). 
So, x(rj +1 ) is the point x(rj) integrated by the flow for a time interval given 

by rj +1 — Tj. Applying the transformation T t j , for an infinite series of rj, that 
is i — 1, 2, 3, . . . , oo, gives us a discrete set that we call V. The projection of 
V on the subspace Vj is named Vj. 

Now, we introduce the notion of transitivity. Let us assume we have a chaotic 
set A. Let us choose two disjoint subsets B and C in A. So, BnC = 0. There is a 
transformation F that generates the set A, with the property that T{B) = C 
and F(C) = B. So, clearly the transformation F and its n-fold application 
(F n ) can always place an arbitrary initial condition, belonging either in B 
or C, anywhere in the set A. This property makes the transformation F to 
be transitive in the set A. However, the conditional Poincare map might no 
longer be transitive. 

For example, in the case of the 2n-fold of the transformation T, i.e., T 2n . Given 
an initial condition in the set B, its iteration by T 2n will never reach the set 
C. Therefore, the conditional Poincare map T 2n is not a transitive transfor- 
mation in the set A. In the case of flows, we have seen that the temporal 
transformation represented by T t j, whenever PS is present, does not place 
arbitrary points of the attractor everywhere in the subsets Xj of the chaotic 
attractor. Therefore, if PS is present the transformation T T i is not transitive 
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to this subset of the attractor. 

The conditional Poincare map T t j is topologically transitive in a set A [2H] if 
for any two open sets B, C C A, 

3rj / T T i\B) nC + 0. (A.l) 



To detect whether the conditional Poincare map is transitive in some subset, 
we introduce the notion of equivalence. Two sets A and B are (not) equivalent 
[equivalence is represented by the symbol =, and non equivalence by ^] if they 
(do not) occupy the same neighboring space. In a more general way: 

Definition 1 Two sets A and B are equivalent, A = B, if Vx G A, a set C 

can be constructed by the union of open sets Bi(x), open R d volumes centered 
at x with length £, such that Wy G B y G C, and A^ B if y ^ C. 

Definition 2 The set T> is a P-set if Xj ^ T>j. 

Definition 3 IfDj can be decomposed into a collection Dj = T>°, T>j , . . . Dj X 
of subsets of Xj, with Q > 1, such that a point in T> 1 - iterated by the conditional 
Poincare map goes to -p* +1 ( mod( 2) ; we re j er f each minimal set as a recur- 
rent decomposition. In the particular case we have rj +1 — rj = r (constant) 
this minimal set is a periodic decomposition. The number of sets Q is called 
length of the decomposition ISUtf . 

Proposition 1 IfT T i is transitive on Xj, then Vj = Xj. 

Proposition 2 If T T j is non transitive in Xj, then a subset A can be con- 
structed, such that A C Xj and A^ Xj. 

Proofs of Propositions 1 and 2 can be done by using Eq. (IA.1I) and the defini- 
tions. 

Definition 3 can be understood by the PS-set of length- 1 and length-2 con- 
sidered in this work. If the length- 1 PS-set is constructed by the time se- 
ries Tj , Tj +1 , Tj +2 , Tj +3 , . . ., the length-2 is obtained by the time series given 

. r i +T i+1 T i+1 +r i+2 

by Tj, — , Tj , — — ^ — ) • • • [25]. For the particular case of a periodically 
forced system, t[ +1 — r\=\jf, and so, the length-1 PS-set is the standard stro- 
boscope map, that collects point every period of the forcing. The length-2 
PS-set is constructed by a stroboscopic map that collects points every half 
period of the forcing. For the length-2 PS-set, we have two minimal sets, the 
sets V°- and V 1 -, with the property that if x G V -, this point iterated by the 
conditional Poincare map goes to V 1 - and V 1 - goes to V -, under the conditional 
Poincare map. 
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To check whether T>j ^ Xj, we do the following, for x G Xj it exists y G T>j 
such that 

ynB e (x)=<&, (A.2) 

where Bi(x) is a open ball of radius 5 centered at the point x. 5 is a small 
positive value. 



B The PS-sets 



In this appendices, we show that a PS-set exists if, and only if, phase syn- 
chronization exists, which in other words PS-set implies PS and PS implies 
PS-sets. 

Given a dynamical system Y' = G(Y), let F* be the flow and X the attractor 
generated by it, we suppose that we have a chaotic dynamics now on. Let 
Ej be the Poincare section in the subspace Vj, and let Uj be the Poincare 
map associated to the section £j, such that given a point x l j G thus 

= rL,'(x*) = F At j {x)). From now on we use a rescaled time if — tj (Ti). 
For a slight abuse of notation we omit the symbol /. 

Proposition 3 Given two interacting oscillators. Then PS-sets can be con- 
structed if, and only if, phase synchronization is present. 

To show the if in the preposition let us start by considering the time interval 
associated to the return of the point to the point x* +1 is Arj +1 = rj +1 — rj, 
with Tj +1 being the times at which the subsystem Xj crosses the Poincare 
section £y, in the subspace Vj. 

Y N _ Ar j 

As already introduced, the average return time is given by (Tj) = — L = 

r N 

4r, and the time is rescaled, such that (Ti) = 1. Our hypothesis is that the 
subsystem Xj has a phase-coherent oscillation, so there is a number 5j for 
which holds [3 lj : 

\T»-N\<8j. (B.l) 

The number 5j < 1 measures the coherence in the phase oscillation, and is 
linked to the phase diffusion [7f31] . This equation holds for all N, so it implies 
that for a single oscillation is also true that |Arj — (T,)| < 5j. In the case of 
two systems that present PS, it holds [21] : 



_N 



T, 



N 



(B.2) 



with £3 < 1. This equation implies \At[ — At^\ < 5%, which states that the time 
intervals in a single oscillation are strongly related in phase synchronization. 
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Let us introduce a new variable that measures the difference between the time 
interval of two events in Xi and X 2 . This new variable is At 2 l = At 2 — Ar{, 
from Eq. fTR2]) . it is true that \At 1 21 \ < S 3 . 

Now, we analyze one typical oscillation. Given the following initial conditions 
%i G Si and x 2 G S 2 , we evolve both until returns to Ei. In other words, we 
evolve both initial conditions for a time Arj. So, F Ar z (x^) = ^(x®) = x\ G 
El Analogously, F A ^(x° 2 ) = F At ' +At ^(x° 2 ) = F Ar l^ oF a ^(x° 2 ) = F At ^{x\). 

Now, we use the fact that (Ar^l < £3, and write that 



F At 2.i (x\) tax\ + G(xl)5 3 . (B.3) 

So, given an initial condition in E x (subsystem X\) evaluated by the time 
initial conditions return in the section E 2 (of subsystem X 2 ), it returns near 
the section E 1; and vice- versa. 

For a general case, we have to show that an initial condition, on the section 
Ei, evolved by the flow for the time Y^=q^ t \\ s ^ m remains close to this 
section. In other words, we have to show that the approximation in Eq. (IB. 31) 
is valid for an arbitrary number of events N in the subspace X 2 . Now, noting 
that X^o^ r 2 1 = T 2 ~ r i 1 fro m our hypotheses of phase coherent dynamics 
I Yh=q ^■ T 2,\\ = \ T ^ N ~ < ^3- The same arguments used to derive Eq. (IB. 3D 
for one oscillation, can be extended to an arbitrary number N, so we proof 
the if of the proposition (PS implies the existence of PS-sets). 

Now, to show the only if (PS-sets imply PS), let us say that there is a PS- 
set. As a consequence, Eq. ( 1B.2j) is valid. Then using the definition of phase 
coherence |31j . and noting that in PS the average return times in a given 
Poincare section are the same, we see that Eq. (1B.2|) comes from the boundness 
of the phase, concluding the proposition. 

Furthermore, let us suppose that the trajectory of the oscillators are perturbed 
by a small perturbation, which does not destroy the phase synchronous dy- 
namics. The effect of the small perturbation in the time return of the trajectory 
to a Poincare section is to deviate this time according Tj = T- + £j, where 
represents the perturbation in system Sj at the moment of the i-th event, with 
maXi\C,j\ < k. Under these hypotheses about the perturbation, we conclude 
the following result. 

Proposition 4 The PS-set is robust under perturbations 

This result shows that a PS-set can be constructed in PS states. Moreover, 
for the coupled Rossler-like systems, this result states that this set is confined 
in an angular region, which is a consequence of Eq. (IB.2H . 
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To see the relation between the constant £3 and the size of the PS-set, we 
do the following. By the time normalization, average time interval between 
points in phase space are proportional to their distances. So, from Eq. (IB. 3D 
we write <5 3 = \Ti\/G{x^), with TL being the average half length of the PS-set. 
A rough calculation shows that in our experiment with the Chua's circuit, we 
have that £3 « 1/2.5, and for the coupled Rossler oscillators, we have that 
#3 pd 1/2.1. This results completely agree with the theoretical approach done 
in [21]. 
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